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One of the outstanding issues in neutrino physics is the experimental determination of the neutrino 
mass hierarchy: Is the order of the neutrino masses "normal" — m\ < m\ < m| — or is it inverted 
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— 7713 < mi < m2, with n%2 — mi <^ m2,mi? While this issue can be resolved in next-generation 
long-baseline ly^ ^ Ue neutrino oscillation studies if |(7e3|'^ = sin^ 613 is large enough, a clear strategy 
on how to resolve it if |(7e3p is sufficiently small is still lacking. 

Here, we study the capability of non-oscillation probes of neutrino masses to determine the neu- 
trino mass ordering. We concentrate on studies of m^^, the kinematical neutrino mass to which 
precise studies of tritium /3-decay are sensitive, ruee, the effective mass to which the rate for neutri- 
noless double-beta decay is sensitive if the neutrinos are Majorana fermions, and E, the sum of the 
neutrino masses, to which cosmological probes of the energy budget of the Universe are sensitive. 
We find that combined measurements of m^e, S, and m^^ are capable of establishing the neutrino 
mass hierarchy if these measurements are precise enough and if one "gets lucky." We quantify the 
previous sentence in detail by performing a numerical analysis of a large number of theoretical data 
sets, for different measured values of m^^, mee, and E, keeping in mind the ultimate sensitivity that 
can be reached by next (and next-to-next) generation experiments. 



I. INTRODUCTION 

In order to explain all neutrino data (with the exception of those from LSND waiting to be confirmed by the 
ongoing MiniBooNE experiment), one is required to augment the standard model Lagrangian by adding operators 
that will render, after electroweak symmetry breaking, the neutrinos massive 

Under the most conservative assumptions, the new standard model parameters that need to be determined from 
experiment are: 

1. The three neutrino masses, nii, i — 1,2, 3, chosen real and nonnegative. These are ordered as follows: ml < m^, 
and < Amf2 = m\ — m\ < jAm^gl, where Amfg = m| — . A positive value of Am^g implies mg > TOj and 
a so-called normal mass hierarchy, while a negative value of Amfg implies rrig < and a so-called inverted 
mass hierarchy. 

The three neutrino masses can be expressed in terms of Amf 2 , Am^g (including its sign) and the lightest neutrino 
mass mi. In the case of a normal mass hierarchy, mi = mi, TO2 — \/ Am^g + ^'^^ '^3 ~ \/ ^"^13 + 
while if the mass hierarchy is inverted m\ ~ — Am^g -I- m^j , m2 = \J — Am^g -I- Am^2 + "'^i^ ^^'^ ™3 — '^i- 

2. The six mixing parameters that characterize the leptonic mixing matrix. We adopt the particle data group 
parameterization of the leptonic mixing matrix (Eq. (13.32) in |j]). Two of the six parameters are unphysical 
if the neutrinos are Dirac fermions. 

Current experiments have been able to measure ^m^^ and jAm^gj, together with the "solar angle" 6*12 and the 
"atmospheric angle" ^23- The third mixing angle 6'ig is constrained to be small (sin^ ^13 < 0.03 at the 99% confidence 
level %\ ) , and we know nothing about the three CP-odd phases 5, ai , 02 . We also don't know whether the neutrino 
mass hierarchy is normal or inverted, and the lightest neutrino mass is only modestly constrained: mi G [0, 2.0] eV, 
at the 99% confidence level Q-* It is among the main goals of fundamental physics experiments to determine the 
currently unknown leptonic parameters. 

In 6], together with Boris Kayser, we explored the issue of determining the mass hierarchy with neutrino oscillation 
experiments. In particular, we discussed the challenges associated with determining the neutrino mass hierarchy in 
the advent that sin^ ^ig is vanishingly small. We concluded that, if ^ig happens to be small enough, the only way 
to determine the neutrino mass hierarchy with neutrino oscillation experiments would be to probe i/^ — > Vy, and/or 
Vy oscillations at very long baselines (L > 5000 km) and relatively small energies [E^ < 500 MeV). Given the 



* This is the most conservative bound on m; , provided by precision measurements of the end point spectrum of tritium beta-decay , and 
ignores stronger constraints from searches for neutrinoless double-beta decay (which apply only if the neutrinos are Majorana fermions) 
and from cosmology (which rely on several nontrivial assumptions regarding the evolution of the Universe and its particle content). We 
will discuss these in more detail later. 
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enormous experimental challenges which one needs to overcome in order to achieve the conditions outlined in 6], it 
is of the utmost importance to explore other probes of the neutrino mass hierarchy. 

Here, we study non-oscillation probes of neutrino masses and mixing, and determine their capabilities for extracting 
the neutrino mass hierarchy. We concentrate on three distinct observables: (i) the electron-type neutrino effective 
mass rriu^ , which can be probed in precise measurements of the end point of nuclear (3 decay spectra, (ii) the effective 
neutrino mass vriee that characterizes the rate of neutrinoless double double beta decay (if the neutrinos are Majorana 
particles), and (iii) the sum of the active neutrino masses S, which is constrained, within the concordance cosmological 
model, by measurements of "cosmological" observables. All three are described in detail in Sec. |n] We spell out the 
circumstances {e.g., how precisely should TOj/^, mee, and E be constrained) under which the observables listed above 
can add to our understanding of the neutrino mass spectrum in Sec. liiil emphasizing the importance of combining 
all three measurements. 

We will concentrate on the limit sin^ 0, when "standard" matter affected <-> v,. searches (see and 

references therein) for the mass hierarchy are known to fail. Under these circumstances, as emphasized above and in 
, alternative probes of the character of the neutrino mass are not only welcome but absolutely necessary. 

II. NON-OSCILLATION PROBES OF THE NEUTRINO MASS 

The only evidence for neutrino masses comes from long-baseline, neutrino oscillation experiments. These take 
advantage of quantum mechanical interference in order to be sensitive to very small neutrino mass-squared differences.* 
On the other hand, neutrino oscillations are impotent when it comes to determining the neutrino masses themselves 
— they are only sensitive to the differences among the various neutrino masses. 

Of course, several other phenomena are also potentially sensitive to nonzero neutrino masses. We discuss three 
qualitatively distinct "non-oscillation probes" of nonzero neutrino masses, and explore their potential for determining 
the nature of the neutrino mass hierarchy under the assumption that it will not be revealed by future oscillation 
experiments. 

A. rrii/^ — Energy Spectrum of Nuclear /3-Decay 

Conservation of energy and momentum dictates that any physical process involving neutrino absorption, emission, 
or scattering is sensitive to the values of the neutrino masses. Due to the smallness of neutrino masses, however, such 
effects are, in practice, almost always hopelessly unobservable. 

Some of the most promising probes are decay processes in which the neutrino energy and momentum are well 
constrained and "easy" to measure with great precision and/or as small as possible. Among these are multi-body 
final state t decays, tt"*" decay at rest, and nuclear beta-decay, which are said to constrain, respectively, the tau-type, 
muon-type, and electron-type neutrino effective masses, which will be properly defined below. 

Given the current constraints on neutrino masses-squared and neutrino mixing, however, the most precise probe, 
by far, of kinematical neutrino mass effects is the precise determination of the end point of the /3-ray energy spectrum 
of tritium decay. The 'end point' is defined as the maximum energy that the daughter /?-ray is kinematically allowed 
to have. 

In a little more detail, the /3-ray spectrum is a function of the three neutrino masses-squared and can be written, 
very schematically, as 

K = \Ue,\'F{ml/ElE,) + \U,2?F{ml/ElE,) -f \U,:,\''F{ml/ElE,), (II.l) 

where K is the /3-ray spectrum as a function of the neutrino energy £;^,1' and F is a function of the neutrino energy 
and the neutrino mass-squared in units of the neutrino energy. In order to obtain Eq. Ijll.l|) . we assume that tritium 
beta decay is described by '^H -^■^ He -f + t'j, and that the probability that a specific mass eigenstate Vi is emitted 

is \U,,\\ 

Eq. Ijll.II) is independently sensitive to all three neutrino masses as long as |C/eiP are known (which is, to a good 
approximation, the case) In the case of very small neutrino masses {mf/E'^ <C 1), however, this is not useful in 



* Hence the need for very long baselines — L 

osc 

It is always amusing to remember that solar neutrino experiments are 

sensitive to Am^ > 10~^^ eV-^. 
t The neutrino energy is trivially related to the energy of the /3-ray, which is experimentally accessible. 
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practice. We can write K, after expanding F around mf/El = 0, as 

K = y: + f ^-i) + o (^) . ("-2) 

+ + (IL3) 

where Fq = F evaluated at m| = 0, Fq = dF/d{mf/El) evaluated at mf = 0, and 

ml=Y,\Ue^\'m^ (11.4) 

i 

is defined as the electron neutrino effective mass-squared. Note that we made use of the fact that J^i \Uei\'^ — 1- 
Henceforth, we will assume, given current bounds on mi and the neutrino oscillation parameters, that tritium beta- 
decay experiments are sensitive only to . 

Experimentally, irif, is currently constrained to be less than (2.0)^ eV^ at the 99% confidence level Given 
that this upper bound is much larger than Ato^j ^'^d jArTi^gl, the current bound on mf is also much larger than 
the neutrino mass-squared differences (this is the so-called "quasi-degenerate neutrino masses" regime) . Under these 
circumstances, ~ ~ rri^ ~ mf, independent of the mass hierarchy, and ~ mf |C/eiP = n^f- Hence, the 
bound TO/ e [0,2.0] eV, quoted in Sec.^ Incidentally, the muon-type and tau-type neutrino effective masses-squared 
are given by to^^ = 'mf\Uai\'^, a = fi,T. Given the current bound to; < 2.0 eV, both are constrained to be less than 
(2.0)^ eV^, orders of magnitude smaller than the current kinematical limits obtained from charged-current processes 
involving muons and taus |3| . 

In the case sin^ ^13 — 0, m^^ is simply given by 

to^ = TO^ cos^ 6'i2 -I- TO.2 sin^ 6'i2, (II-5) 

= TO^ + Am^2sin^6'i2, (II. 6) 

— TO^ -I- Ato^2 sin^ 012 (normal hierarchy), (II-7) 

= mf — Ato^3 -I- Am^2 sin^ ^12 (inverted hierarchy). (II-8) 

In the case of a normal mass hierarchy, the expressions above still hold for all values of to/ and nonzero \Ue3\, as long 
as I Ato^3| sin^ 6*13 <^ Ato^2 sin^ 6*12 — > sin^ ^13 <C 0.01. In the case of an inverted hierarchy, the expressions above 
are valid if sin^ 6113 <C 1, which is already experimentally guaranteed to be satisfied. Fig.^depicts m^^ as a function 
of m; for an inverted and a normal mass hierarchy, and for Ato^2 ~ 8.0 x 10~^ eV^, Ato^3 = —2.44 x 10"'^ eV^, 
sin'^ 012 — 0.31, and (as already mentioned) sin'^ ^13 — 0. 

We are interested in extracting the neutrino mass hierarchy from experimental information on to^^, and address 
the question in the following way. We assume that all oscillation parameters are measured precisely enough that 
uncertainties on these are irrelevant (we worry about the precision with which oscillation parameters are measured 
in the next section). Furthermore, since we are working under the assumption that |C/e3| = 0, we also assume that 
the neutrino mass hierarchy is not known, i.e., there are two values of Ato^3 that perfectly fit all neutrino oscillation 
data. These will be referred to as Ato^;^ > and Ato^3 < 0. As discussed in detail in 0], it is quite likely that 
Ato^J^ — |Ato^3 I = X ^ 0. The value of x will depend on the details of the measurement of Ato^3 

Currently, the experimental upper bound on m^^ does not allow one to discriminate a normal from an inverted 
mass hierarchy. This remains true as long as 

ml^ > -Ato^3 . (II.9) 

Indeed, given our current knowledge of Am^^, in order to determine the mass hierarchy from an upper bound (or 
measurement) of the electron-type effective neutrino mass, one should be able to constrain < 1.7 x 10^"^ eV^ (the 
current 99% upper bound on Am^;^ 5]). This can only happen if the mass hierarchy is normal. 

In other words, if the neutrino mass hierarchy is inverted, it is not possible to determine the character of the 
neutrino mass by measuring m^^ , irrespective of how precise a measurement one is able to perform. Associated to 
the inverted mass hierarchy measurement of mi = mf , there is a different value of to/ = m^ that provides an equally 
good fit to all neutrino data as long as mf^ — mf~ — Ato^;^ (we remind readers that Ato^^ is negative-definite, 
such that m^ > mf). If the mass hierarchy is normal and the true value of the lightest mass is to^^ > — Ato^^^, 
there is an equally good fit to all the neutrino data with an inverted hierarchy and a lightest mass mj given by 
TO; ~ = mi + Atoj^3 . If TO; < — Ato]^!^, there is no inverted mass hierarchy "mirror solution." 
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FIG. 1: m^^ (dotted lines) , 
mass hierarchy (normal A'^ or inverted /) for Ami2 - 
sin^ 9i2 = 0.31, and sin^ ^13 = 0. See text for details. 



(solid bands) and E (solid lines) as a function of the lightest neutrino mass, mi and the neutrino 



8.0 X 10"^ eV^ Am?+ = 2.50 x 10"^ eV^ Am? 



-2.44 X 10"^ eV^ 



In the future, the KATRIN experiment, currently under construction, aims at being sensitive to values of m^^ > 
0.2 eV at the 90% confidence level Q. Note that this stiU satisfies Eq. ||IL9|) . such that KATRIN will not be able 
to establish the neutrino mass hiearchy, even if it observes a nonzero m^^ effect. We are not aware of experimental 
proposals to significantly improve the sensitivity to nii,^ beyond the reach of KATRIN. 



B. TUee — Neutrinoless Double Beta Decay 



If the neutrinos are Majorana fermions, lepton number is not a good quantum number. -I- The most powerful probe 
of lepton number violation is the search for neutrinoless double beta decay, Ovfifi: [Z + 2)e^e^. Assuming 

that the neutrino masses are the only relevant source of lepton number non-conservation, the rate for Qvj3[3 is governed 
by the fundamental process W~W~ — > e~e~, mediated by Majorana neutrino exchange. 

Schematically, 



rrii 



(11.10) 



where is a typical energy momentum transfer (squared) associated with the OvPP process — ^ (50)^ MeV^ — 
and M. contains the rest of the kinematics and the nuclear matrix element associated with the Z — > (Z + 2) transition. 
As is well known, the rate for Oi//3/3 is proportional to a combination of the neutrino masses that vanishes in the limit 
rrii —>■ 0, Vi. This is, of course, expected, since we are assuming that nonzero Majorana neutrino masses are the only 
source of lepton number violation. 



t To be more precise, U{1)b—l is no longer a global symmetry of the Lagrangian that describes standard model degrees of freedom. 
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+ ;^)), (11.11) 



All m| are much smaller than Q^, so we can simplify Eq. (jlLlOp by expanding the neutrino propagator 
and M{mf, Q'^) — M{0, Q^){1 + 0{mf /Q^)) such that, up to corrections proportional to mf/Q"^, 

roj./3/3 OC |meeP, (11.12) 

where 

mee = ^(7^™, -milL/eipe^"^ +TO2|t/e2pe*"' +TO3|C/e3pe-2^^ (11.13) 

i 

is the effective mass for 0^(3(3.^ We will assume henceforth that decay experiments are capable of constraining 

(or measuring) |mee|. 
In the case of interest here (sin^ (?i3 — 0), 



and, 



|mee| = |miCos2 6'i2e'"i +m2sin2 6li2e'"^| (11.14) 

1 /2 

= (rn^ cos^ 6'i2 + 7712 sin^ ^?i2 + 2TOim2 sin^ 012 cos^ 012 cosa) ' (11.15) 

/ sin^ 29 \ 

mf cos"* 012 + (mf + Ami2) sin'' 0i2 + miJmf + Amfj cos a I (normal hierarchy), (II. 16) 



{[nif - Arn^g) cos'* 0i2 + [mf + Ato^2 - Am^g) sin'' 0i2 

-sin2 2012 ^ ^'^ 



^ (mf — Aml^){mf + Ami2 — Amig) ^^^ — ^cosaj (inverted hierarchy), (11-17) 



where a = a2 — ai is the Majorana phase to which rriee is sensitive (in the Hmit sin^ 0i3 — 0). Given the current 
constraints on sin^0i3, the expressions above are an excellent approximation for niee in the ca se of a n inverted 
hierarchy, and also apply safely for a normal mass hierarchy as long as a/| Am^gj sin^ 0i3 ^ y^Am^sin^ 0i2 
sin^ 013 ^ 0.05. Fig. ^ depicts \mee\ as a function of m; for an inverted and a normal mass hierarchy, and for 
Am^2 = 8.0 X 10"^ eV^, Am?3 = -2.44 x lO'^ eV^, sin^ 0i2 = 0.31, and sin^ 0i3 = 0. Associated to a fixed value of 
mi and a mass hierarchy there is a continuum of values of \mee\, one for each value of cos a € [—1, 1]. 

The minimum (maximal) value of ruee for a fixed value of to/ corresponds to cos a = —1 (+1). rriee can vanish in 
the case of a normal mass hierarchy, when mf cos'' 0i2 = (Ami2 + mf) sin* 0i2. Hence, for 

TO? ^ Ami, ,f"'^"\, - (0•004)^ (11.18) 

COS^ 012 — Sm 012 



and an inverted mass hierarchy, m^e is especially suppressed. 
As in the previous subsection, it is easy to show that if 



iTOeel > \l Aml^ cos^ 0i2 - ^ ATO^g + Ami, ^ 0-02 eV, (11.19) 

a measurement of |TOee| has no discriminatory power when it comes to determining the nature of the neutrino mass 
hierarchy, irrespective of its precision. As before, such values can only be experimentally ruled out if the mass hierarchy 
is normal ■9|. Hence, any "perfect" fit to all the neutrino data (including a potential measurement of Iweel) obtained 
under the assumption that the neutrino mass hierarchy is inverted can be matched by an as-perfect fit obtained 
under the assumption that the mass hierarchy is normal. A very recent and thorough analysis of the capabilities of a 
measurement of \mee\ to determine the mass hierarchy can be found in [lOj . 



^ After electroweak symmetry breaking and in the weak basis where the charged current couplings and the charged lepton mass matrix 
are diagonal, rriee is the 'Vei^e" entry of the neutrino Majorana mass matrix. 
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In practice, the situation is significantly more involved, for three main reasons. One is that the extraction of |mee| 
from the rate for Oi//3/3 is severely clouded by uncertainties in computing the nuclear matrix elements. Different 
theoretical estimates that make use of different nuclear physics techniques can differ significantly, often by an order 
of magnitude HjJ. We have nothing to add to the discussion, except for the fact that the situation should improve 
significantly in the next several years, especially if a positive signal for QvPP is obtained. Another source of confusion 
is that Tq^pp is only proportional to \mee \ if there are no other "beyond the standard model" sources of lepton number 
violation. We will come back to this issue briefly in the next section, and argue that the combined analysis of searches 
for rriv^, jweel and S (discussed in the next subsection) can shine light on this subject. Finally, the rate for Qvj3[3 is 
only nonzero if lepton number conservation is not exact. As far as this discussion is concerned, this implies that while 
evidence for Qv(3(3 can be translated into a nonzero |mee|, failed attempts to observe Qv(5j3 cannot be translated into 
an upper bound on |mee|, and hence potentially vital information regarding the nature of the neutrino mass hierarchy. 
This should be contrasted with failed attempts to measure rrii,^ from tritium beta decay. These can be (much more 
model-independently) translated into information that may help resolve the nature of the neutrino mass hierarchy. 

Currently, niee is (conservatively) constrained to be |mee| < 0.91 eV at the 99% confidence level (see and 
references therein for details) , and we will ignore the to-be-confirmed recent evidence for a nonzero rate for OvjSfi |l2j . 
As in the case for the current bound on mi^^ , no information regarding the neutrino mass hierarchy can be obtained 
from our current knowledge of |TOee|- This would remain true even if we postulated that the neutrinos are Majorana 
fermions. Near-future experiments with sensitivity to |mee| 0. 1 eV are currently being planned, and several of those 
aim at reaching an ultimate sensitivity to |mee| ^ 0.01 eV It is fair to say that significant R&D efforts (plus 
time and resources) are necessary in order to improve the sensitivity to |mee| beyond that. 



C. E — Cosmological Observables 



There is a significant amount of indirect evidence that the Universe is filled with a very cold thermal population 
of relic neutrinos — the cosmic neutrino background. While we have no direct evidence that these really exist, their 
presence is required in order to match predictions of big bang nucleosynthesis with measurements of the primordial 
abundances of light nuclei Indeed, before LEP experiments measured the invisible .Z-boson width, bounds from 
big bang nucleosynthesis were known to provide the most stringent constraints on the number of light neutrino species 



If neutrinos have mass, the cosmic neutrino background contributes to the energy budget of the Universe as hot 
dark matter. According to the very successful (if not very palatable) concordance cosmological model 15], the fraction 
of hot dark matter in the Universe is constrained to be small. Given our current understanding of neutrino masses, 
rihot — X] "^i/(48 eV), such that an upper bound on fihot can be translated directly into an upper bound on S = m^. 
As was the case in previous subsections, cosmological observables are potentially sensitive to more than the sum of 
the neutrino masses, at least in principle. For example, hot dark matter has the property of preventing matter from 
"clumping" in the early Universe, leading to a suppression of power at small scales in the matter distribution power 
spectrum [T^ . Potential "kinks" in the spectrum are related to the mass of the individual hot dark matter species. 
We will ignore this possibility and will assume that cosmological observables are only sensitive to S. 

S depends only on the neutrino masses, and can be expressed in terms of m; and the mass-squared differences: 



S = mi + i/mf Am|2 + Wmf + Amfg, (11.20) 



m; -t- w 771^ + Amj2 + y"*; + Am^3 (normal hierarchy), (11-21) 



= mi + \J'mf — Am^3 + y mf + Am^2 ~ ^"^13 (inverted hierarchy). (11.22) 

E is the only observable discussed in this section that, in the limit ^13 0, is sensitive to Am'l'^ , the value of the 
atmospheric mass hierarchy obtained assuming that the mass hierarchy is normal. Fig. ^ depicts E as a function of 
the lightest neutrino mass, mi for Amfa = 8.0 x 10"^ eV^, and Aml^ = 2.50 x 10"^ eV^ Aml^ = -2.44 x 10"^ eV^ 
for both neutrino mass hierarchies. 
As before, if 



E > ^J-Aml^ + y Amfa - Am^^ , (11.23) 

a measurement of E is not sensitive to the nature of the neutrino mass hierarchy, independent of its precision. In 
light of the discussion in the two previous subsections, it should not be surprising that Eq. Ijll.23p is automatically 
satisfied in the case of an inverted mass hierarchy, but could be violated in the case of a normal one, as long as mi 
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is small enough (m/ ^ 0.02 eV). Quantitatively, for every candidate value of the lightest mass mi = mj' obtained by 
postulating that the mass hierarchy is inverted, there is another lightest mass value mi = that fits the neutrino 
data just as well as long as one postulates that the mass hierarchy is normal and satisfies 



ml + y mf+ + Am^2 + V "'-f^ + ^™13 = ™; + V - ^"^13 + V "^f + ^^-12 ~ ^"^13 ■ (11.24) 

As with measurements of |mee|, the map between cosmological measurements and E is nontrivial and model 
dependent. The cosmological determination of S depends on the history of the Universe and its particle content, 
and physics beyond the standard model and/or concordance cosmology can easily obscure, enhance or completely 
erase any relationship between cosmological observables and E 0, 0| . 

According to j^, combined analysis of data from WMAP, SDSS and Lyman-alpha forest surveys constrain E < 
0.94 eV at the 99% confidence level. Different analyses of the same data or subsets thereof quote similar bounds 
(within 50%) 1^. Upcoming missions have the potential to significantly improve on the current sensitivity to a 
nonzero value of E. It is expected that future data on the cosmic microwave background should be sensitive to 
E > 0.1 eV while current studies of the capabilities of weak lensing probes seem to indicate that these might 
be sensitive to E values below the 0.1 eV level For an overview see, for example, [l^ . It has recently been 

speculated that one may indeed reach E values as low as 0.03 eV (201] . 



III. COMBINED FITS TO THE MASS HIERARCHY 

It has been widely recognized in the literature l2lll2^ I23L l2^ that the combined analysis of searches for rriee,* 
and/or E provide qualitatively more inforjnation than the separate analysis of each observable, including potential 
information regarding the mass hierarchy |2^ . For example, a lower bound on m^,^ combined with an upper bound for 
may help determine, with the addition of some reasonable assumptions, that the neutrinos are Dirac fermions, 
while a large m^^ value, combined with an upper bound on E would lead one to conclude that there is more to 
our understanding of particle physics or the history of the Universe than currently accepted. Here, we discuss the 
circumstances under which measurements of (or upper bounds for) m^^, m^e, and/or E help us determine the mass 
hierarchy. 

It is easy to understand why a combined analysis should prove more powerful. The value of that renders a normal 
hierarchy fit to, say E, compatible with an inverted hierarchy one (associated to a value mf of the lightest mass) need 
not agree with the the equivalent quantity obtained from analyzing measurements of m^e or m^^ . For example, if the 
true value of mi vanishes and the mass hierarchy is inverted, the mirror value mi = to^ that renders E "the same" if 

one assumes that the mass hierarchy is normal is m^(E) ^ O.Sy^Aro^ ~ 0.02 cV. On the other hand, the equivalent 

for m^^ is m^i^m^^) — \J ^m\^ ~ 0.05 eV. Hence, for the wrong-hierarchy hypothesis, mi{m^^) ^ mi(E) ^ mi{mee) 
— even if the true neutrino mass hierarchy is inverted. This indicates that combined analyses of the three observables 
discussed in Sec. ^] can improve, sometimes significantly, our ability to determine the neutrino mass ordering. 

Figure 121 depicts the effective masses pairwise for both hierarchies, assuming Am^2 = 8.0 x 10~^ eV^, Am^j^ = 
2.50xlO~^eV2, Amlj = -2.44x 10"^ eV^, sin^ 6*12 = 0.31 and sin^ 6*13 = 0, allowing rri; to vary between and 0.5 eV, 
and the relevant Majorana phase to span its entire allowed physical range: cos a G [—1, 1]. One can easily identify 
sets of potential measurements that, with sufficiently small uncertainties, would imply a specific mass hierarchy.^ In 
particular, it is easy to verify, and improve on, the claims made in section ^] regarding the implication of the normal 
scheme for sufficiently small measurements of the effective masses. Evidence for the inverted hierarchy can also arise 
from such naive considerations. For instance, by examining the upper left panel it is clear that the inverted hierarchy 
could be established if accurate measurements of rriee and E were made near 0.04 eV and 0.1 eV, respectively. The 
same conclusion is obtained if only an upper bound exists for E near 0.2 eV, while a measurement sets mee near 
0.05 eV. It is easy to see that combined measurements of m^^ and mee cannot distinguish an inverted mass hierarchy 
from a normal, quasi-degenerate scheme — knowledge of E is extremely valuable when it comes to establishing an 
inverted mass hierarchy. 

This simple picture is smeared out when realistic conditions are taken into account. The inclusion of parameter 
uncertainties will naturally wash out some of the interesting regions in this measurement space, but most of the main 



* Henceforth, to simplify the notation, we define rriee = |mee|. 

t The usefulness of looking at the non-oscillation observables as depicted in Fig. |2] was pointed out in .:23l . 
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FIG. 2: Values of m^^, |mee|, and E for both hierarchies, for Ami2 = 8.0 x 10 ^ eV^ 
-2.44 X 10"^ eV^ sin^ 6I12 = 0.31 and sin^ 6113 = 0. 
phase spans its entire allowed physical range 



2.50 X 10"^ eV^ 



The value of mi varies between and 0.5 eV and the relevant Majorana 
cosQ £ [—1, 11. See text for details. 



features remain robust. A complete statistical analysis, quantifying the very visual and intuitive line of reasoning 
outlined here, is presented below. 



A. Parameters and Data Analysis 



In order to quantitatively compare the different neutrino mass hierarchies we proceed as follows. We assume that all 
oscillation parameters Qi = Am^s , Am^2 , sin^ 6\2, sir? 613 are measured to be 8i ± ae; , while the three non-oscillation 
parameters are constrained to be m^^ = fn^ ± a^^ , rUee = Wee i Cee , and S = E ± cts . Note that we are assuming 
that the uncertainty related to the computation of the Qv(3[3 nuclear matrix element is included in Gee ■ We assume 
no information concerning the three CP-odd phases ai,a2,(5. We refer to all these measurements as the data. We 
further postulate that the neutrinos are Majorana fermions, and that there are no other sources of lepton number 
violation. Finally, we assume concordance cosmology and "standard" particle physics. 

Armed with these results, we construct a function defined by 



X^(m;,ai,a2,(5, 9) = 



[rriy^ - myXmuQ)Y' [rTiee - TOee(TOi, 6, ai, Q!2, (5)]^ [S - E(r7i;, 9)]^ 



E 



(ni.1) 

This function also depends implicitly on the choice of hierarchy from the structure of the effective masses m^^, mee, 
and E given, respectively, by Eqs. HII.4III.13III.20p . 

For each hypothesis concerning the mass hierarchy, we minimize with respect to all neutrino parameters {Qi, 
mi, and the two relevant CP-odd phases [say, a2 — ai, and —25 — ai]), and establish if a good fit to the data can 
be obtained. This is done by comparing the minimum value of x^ = Xmin with expectations for three degrees of 
freedom. In particular, Xmin < '^■^ implies that the theory fits the data at the 95% confidence level. If only one of the 
hypotheses concerning the mass hierarchy fits the data at the 95% confidence level (or better), we conclude that it is 
established (at least) at the 95%, while if both hypotheses prove to fit the data at the 95% confidence, we conclude 
that the mass-hierarchy cannot be determined. Finally, if neither hypothesis fits the data, we are forced to conclude 
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TABLE I: Input values, including uncertaintites, for the simulated data sets, 'x-axis' and 'y-axis' refer to Fig.|3 The 'comments' 
column refers to the experiment (or class of experiments) expected to contribute most significantly to the uncertainty with 
which each neutrino observable is assumed to be measured. 



Parameter 


Run 1 
Central a 


Run 2 
Central a 


Run 3 

Central g 


comments 


Am?2(10"^ eV^) 
Am?+(10-^ eV^) 
Amis (10~ eV^) 

sin^ 6^12 

sin^ 6i3 


8.0 5% 
2.5 1% 
-2.44 1% 
0.31 5% 
0.01 


8.0 5% 
2.5 1% 
-2.44 1% 
0.31 5% 
0.01 


8.0 5% 
2.5 1% 
2.5 1% 
0.31 5% 
0.01 


Reactor [29] 
T2K, ^Ovk, MINOS [2^ 
T2K, NO;/A, MINOS [2^ 
Reactor [28] [29] 
T2K, NO;/A,D-CHOOZ [25] 


m,j^ (eV) 
rriee (eV) 
E (eV) 


0.1 
y-axis 0.01 
X-axis 0.01 


0.1 
y-axis 0.05 
X-axis 0.05 


0.1 

y-axis 0.1 
X-axis 0.1 


KATRIN |8] 

m 







that there is a flaw in our theoretical model. Candidate flaws include the hypothesis that the neutrino is a Majorana 
fermion, that concordance cosmology is correct, etc. 

The different data sets used in our numerical analysis are listed in Tabled The best fit values Q from the global 
analysis of neutrino oscillations were used for 0, except for Am^g. In the case of a normal mass hierarchy, we pick 
Am^+ = 2.50 X 10^3 gy^, while in the case of a normal mass hierarchy, we choose Amjg = — Am^^ -I- where 
X = 2Ato^2 cos 2^12, in order to take into account that fact that, in the case of very precise measurements of Amfg, 
different hypothesis regarding the mass hierarchy yield different values for jArn^g] The value of x above agrees with 
what one expects to obtain in accelerator studies at relatively short baselines {L < 3000 km) or relatively 

large energies (Ei^ > 1 GeV). We refer readers to Q for details. 

In order to anticipate future developments, we assume that the uncertainties on the various oscillation parameters 
are smaller than the current ones. For instance, the uncertainty on the atmospheric mass-squared difference ]Amf3] 
should be reduced to ~ 3% in the future, due to potential results from off-axis beam experiments, such as T2K and 
NOi^A 1251 ]. and future endeavors can improve on that 26]. Similarly, we expect that the the solar parameters, A?7i|^^ 
and sin 6*12 can be measured at the 5% level in, say, precise measurements of the low energy solar neutrino flux [23, 
improved reactor experiments boasting greater statistics and reduced systematic errors conducted at baselines at (or 
near) the first survival probability minimum |2^ |2^ , or with broadband long-baseline v^. studies ,2g,] . 

B. Results 

Figure|31depicts regions of the S x plane consistent with both a normal and an inverted neutrino mass hierarchy 
(green [gray]), consistent only with a normal mass hierarchy (yellow [lightest gray]), consistent only with an inverted 
mass hierarchy (blue [dark]), or inconsistent with either hypothesis concerning the mass hierarchy (white). The 
uncertainties with which S and are expected to be measured, o^^^ and cts, are held fixed for all values of E, 
rn^, and are both equal to 0.01 eV in "Run 1" (left panel), 0.05 eV in "Run 2" (middle panel) and 0.1 eV in "Run 3" 
(right panel). Throughout, we assume that is experimentally constrained to be less than 0.1 eV (at the one sigma 
level).''' In the absence of new sources of lepton number violation, this upper bound on implies Jmee] ^ 0.1 eV 
and, in the absence of nonstandard cosmology, S % 0.3 eV. We concentrate our exploration of the the parameter 
space to this region. If, on the other hand, the KATRIN experiment obtained statistically significant evidence that 
m^^ is nonzero (say, m^^ = 0.3 ±0.1 eV) we would be constrained to the quasi-degenerate neutrino mass scenario, in 
which it is virtually impossible to determine the mass hierarchy by precisely measuring m^^ , m^^, , and S (see Fig.O)). 

If both S and are only poorly measured (cTmee ='^11 = 0.1 eV, Fig. EIRIGHT)) it is, in general, not possible 
to determine the mass hierarchy, irrespective of S and mee- The only exception corresponds to very large values of 
mee, where, given an upper bound on E, only an inverted mass hierarchy provides a good fit. Very large values of E 



■f This is in rough agreement with the sensitivity of the KATRIN experiment, which can exclude < 0.2 eV at the 90% confidence 

level if it does not observe any evidence for nonzero neutrino masses 
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FIG. 3: Regions of the E x rriee plane consistent with both a normal and an inverted neutrino mass hierarchy (green [gray]), con- 
sistent only with a normal mass hierarchy (yellow [lightest gray]), consistent only with an inverted mass hierarchy (blue [dark]), 
or inconsistent with either hypothesis concerning the mass hierarchy (white) . Input parameters and expected uncertainties are 
tabulated in Table 13 LEFT: cr„„ = cte = 0.01 eV, MIDDLE: cr„„ = (te = 0.05 eV, and RIGHT: a™^^ = cte = 0.1 eV. See 
text for details. 

and are inconsistent with the upper bound on m^^ , and, if observed, need to be interpreted as evidence for new 
cosmology/particle physics and sources of lepton number nonconservation other than the neutrino masses. If E and 
iriee are measured somewhat more precisely (cr„j_._. = (t^, = 0.05 eV, Fig. I^JMIDDLE)), the situation is qualitatively 
the same: one will not be able to determine the mass hierarchy from non-oscillation probes of neutrino masses unless 
niee is measured to be large (m^ 0.1 eV) and E is "small" (E < 0.1 eV), in which case only the inverted mass 
hierarchy provides a good fit to the data. 

Only when both S and TOgg are measured very precisely {u^^^ = o^s = 0.01 eV, Fig. |3|[LEFT)) do we start to 
also positively discriminate a normal mass hierarchy. If E ;< 0.07 eV, only the normal mass hierarchy hypothesis is 
capable of fitting the data, and only if rriee ^ 0.04 eV. As before, large m^e values, appropriately correlated with 
large E values, point to an inverted mass hierarchy. Furthermore, when the sensitivity of both probes reaches the 
0.01 eV level, we fail to fit the data if Wiee ^ 0.03 eV and E > 0.3 eV. If faced with such a scenario, we would be 
tempted to conclude that the neutrinos are Dirac fermions. Finally, note that, if as — 0.01 eV, we are guaranteed, in 
the absence of new cosmology /particle phisics, to obtain a nonzero value for E — hence the region in Fig. OlRlGIIT) 
corresponding to E < 0.02 eV, for all m^, is white. 

As alluded to above, if ruee ^ 0.1 eV and S ^ 0.3, we enter the quasi-degenerate neutrino mass regime, where 
discrimination between the two mass orderings is virtually unachievable. This behavior can be clearly observed in all 
panels of Fig. 13 

IV. CONCLUDING REMARKS 

While we have discovered that neutrinos have mass, there are still qualitative aspects of the neutrino mass spectrum 
that remain unknown. We have been able to measure, with good precision, the (absolute value of the) two mass- 
squared differences, but know very little about the magnitude of the individual neutrino masses. We are also ignorant 
when it comes to the neutrino mass hierarchy. We don't know if the neutrino masses are "normal ordered" — 
< m| < m§ — or whether the mass hierarchy is "inverted" — m§ < m\ < rn^ (such that Amfj m\, m\). The 
reason we have failed (so far) to discover the neutrino mass hierarchy via neutrino oscillations is that both [J/esI^ and 
Am^2/I ^"^isl happen to be small. 

Next-generation, long-baseline neutrino oscillation experiments are poised to uncover the ordering of the neutrino 
masses, but they can only be successful if |?7e3p is large enough. If, on the other hand, |f/e3p is too small, it is quite 
challenging to determine the neutrino mass hierarchy via neutrino oscillations 0. Other probes are necessary, kere 
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we have studied in detail whether combined information from future studies of neutrinoless double beta decay (mge), 
tritium beta-decay (to^^) and cosmological probes of the energy composition of the Universe (E) can help establish 
the neutrino mass hierarchy. 

Our results are summarized in Fig. |21 and, as expected, the situation is quite challenging. Assuming we can rule 
out (at the one sigma level) rriu^ > 0.1 eV and probe TOge and S with Um^^ 0.05 eV, we can only hope to 

determine the neutrino mass hierarchy from nonoscillation experiments if it is inverted and if the Majorana phase to 
which TTiee is sensitive is close to ("constructive interference," cosa ~ 1).* If cTmee ~ crs ~ 0.01 eV, the situation 
is significantly improved, and there is also the possibility of establishing a normal mass hierarchy if S < 0.07 eV 
and TOee ^ 0.04 eV. In the next several years (^ 10), it is reasonable to expect that we will be somewhere between 
Fig-ElMGHT) and Fig. |3J MIDDLE), while it is not overly optimistic to assume we will eventually reach somewhere 
between Fig. |3;MIDDLE) and Fig. |3;LEFT). 

It is curious to note that the "main goal" of non-oscillation neutrino experiments is, arguably, not to determine the 
mass hierarchy. These experiments are usually associated with establishing if the neutrinos are Majorana fermions, 
measuring the lightest neutrino mass mi, and, ultimately, probing whether the Majorana phases are nontrivial 
(ai,a2 7^ 0,7r). While we have not concentrated on these issues here, it is clear that in the green [gray] regions 
of Fig. 121 one is not able to properly measure m;, even when niee and E are measured to be nonzero with good 
precision. In those regions, for every normal-hierarchy measurement of nii = ni^ ± , there is an inverted hierarchy 
mi = mj' ± that fits the data just as well. 

We reemphasize that, in order to convert measurements of Foiz/j/j and cosmological data into information regarding 
neutrino masses, we are required to rely on several untested hypothesis, including the fact that Majorana active 
neutrino masses are the only source of Icpton number violation and the fact that we quantitatively understand the 
history of the Universe from times slightly before big-bang nucleosynthesis until the formation of the observed large- 
scale structure. It is reassuring that combined measurements of m-j,^, mee, and S can help falsify these hypotheses in 
some regions of the parameter space, but there is a finite possibility that, after many years of data analysis, we will be 
left in a rather "confused" state. This would happen, for example, if the mass hierarchy were normal, mi < 0.01 eV 
and (Ts ^ 0.03 eV. 

We conclude by pointing out that, as far as establishing the neutrino mass ordering is concerned, we would profit 
tremendously with more precise information on m^^ . A measurement of (or upper bound for) m,y^ at the sev- 
eral xlO~^ eV level, combined with a^, slightly below the 0.1 eV level would be sensitive to the mass hierarchy (see 
Fig- ! and it is easy to appreciate that extracted values of m^^ are significantly less model dependent than those of 
mee and S. 

Note Added: While this work was being completed, was posted in the preprint ArXiv's. It shares several of 
the results presented here (with a slightly different treatment of the "data" ) , and contains a very detailed discussion of 
the capabilities of T, and mee (including a detailed treatment of the uncertainties related to nuclear matrix elements) 
to establish the neutrino mass hierarchy and the existence of nontrivial Majorana phases, and to measure to/. Our 
contribution, on the other hand, contains a more detailed analysis (including a detailed qualitative discussion) of the 
capabilities of combined m^^, mee, and S measurements to uncover the neutrino mass ordering, motivated by the fact 
that it will remain unknown if | Ue3 \ is vanishingly small, and also discusses the importance of measuring (or further 
constraining) m^^. 
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